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ABSTRACT 

We present infrared bolometric luminosity corrections derived from the detailed spec- 
tral energy distributions of 62 bright quasars of low- to moderate-redshift (z = 
0.03 — 1.4). At 1.5, 2, 3, 7, 12, 15, and 24 /mi we provide bolometric corrections 
of the mathematical forms Li SO = £ \L\ and log(Lj so ) = A + B \og(\L\). Bolometric 
corrections for radio-loud and radio-quiet objects are consistent within 95% confi- 
dence intervals, so we do not separate them. Bolometric luminosities estimated using 
these corrections are typically smaller than those derived from some commonly used 
in the literature. We investigate the possibility of a luminosity dependent bolometric 
correction and find that, while the data are consistent with such a correction, the 
dispersion is too large and the luminosity range too small to warrant such a detailed 
interpretation. Bolometric corrections at 1.5 /xm are appropriate for objects with prop- 
erties that fall in the range log(L{, ;) = 45.4 — 47.3 and bolometric corrections at all 
other wavelengths are appropriate for objects with properties that fall in the range 
log(L 6o/ ) = 45.1 - 47.0. 
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1 INTRODUCTION 



X 



Recently in iRunnoe et all |2012l ), we derived optical, ultra- 
violet (UV), and X-ray bolometric luminosity corrections 
fro m the quasar spect ral energy distributions (SEDs) of 
the IShang et all (|201 lh atlas. We addressed the practical 
difficulties of determining bolometric corrections, including 
source variability, contamination by emission from the host 
galaxy, gaps in data coverage, and the subtleties of assuming 
isotropy when integrating the SED, and derived bolometric 
corrections at 1450, 3000, 5100 A, and 2-10 keV. The 
bolometric corrections we provided are available with zero 
and nonzero intercepts, are appropriate for radio-loud (RL) 
and radio-quiet (RQ) quasars, and come with an additional 
correction for viewing angle that is not built in. Accretion 
discs do not emit isotropically; we differentiated between 
the bolometric luminosity calculated under the assumption 
of isotropy and the true bolometric luminosity which mo- 
tivated the additional correction to be applied after their 
optical/UV or X-ray bolometric corrections. 

In April of 2011, the first public data rel ease from the 
Wide -field Infrared Survey Explorer (WISE; IWright et all 
|2010| ) opened wide a window into the infrared (IR) of active 
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galactic nuclei (AGN) that was most recently explored by 
the Spitzer Space Telescope (I Werner et al.ll2004l ). The popu- 
lation of IR-selected AGN that is revealed by these missions, 
often reddened quasars, may be unavailable at optical wave- 
lengths due to obscuration and the refore cannot u s e the op- 
tical/UV bolometric corrections of IRunnoe et all (|2012t ) as 
a viable method for determining bolometric luminosity. 

There are infrared bo l ometr ic corrections available in 
the literature. lElvis et all (| 19941 ). the original SED atlas, 
provides a bolometric correction at 1.5 /im, b ut the SEDs 
are no t as detailed and up-to-da te as those of Shang et all 
|201ll ). especially in the infrared. iRichards et al-H 20061 ) em- 
ploys a much larger sample size in exchange for less complete 
data coverage and provides a bolometric correction at 3 fim. 
Both of these corrections ar e subject to the effec ts of dou- 
ble counting in the infrared. iHopkins et ail l|2007fl derives a 
double po wer law correction at 1 5 fim that is broadly consis- 
tent with IRichards et all |2006l ) from an observation-based 
model SED. 

With the release of WISE and the ensuing level of in- 
terest in the infrared emission of quasars, there is a need for 
new infrared bolometric corrections. W e determine infrared 
bolometric corrections from the S EDs oflShang et all d201ll ) 
using the same methodology as in IRunnoe et all (|2012t ). 

This paper is organized as follows. Section 2 describes 
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the data and the lShang et all (|201lh SEDs. In Section 3 we 
derive bolometric corrections and in Section 4 we discuss 
them in the context of previous work and recommend the 
best bolometric corrections for various situations. Section 5 
summarizes this investigation. 

We adopt a cosmology with Ho = 70 km s _1 Mpc -1 , 
SIa = 0.7, and Q m — 0.3. Note that we also differentiate 
between the bolometric luminosity calculated under the as- 
sumption of isotropy (L iao ) and true bolometric luminosity 
(Lboi), which likely differ. 



2 SAMPLE, DATA, AND MEASUREMENTS 

Our empirical bolometric corrections are deri ved from the 
infrar ed-to-X-ray continua of a subset of the IShang et al.l 
l)201ll l SEDs. We give a brief summary of the SED sam- 
ple , data, and bolome tric l uminosities, wh i ch are discussed 
by IShang et all l^OUl ) and iRunnoe et"al] i|2012l) in detail. 
We then present measurements of monochromatic infrared 
luminosities from the SEDs. 



2.1 The SED sample 

The atlas has a total of 85 objects from three different sub- 
samples. The 'PGX' subsample consists of 22 of 23 Palomar- 
Green (PG) quasars i n the complete sample selected by 
lLaor et al] (|1994l . ll997h to study the soft-X-ray regime. This 
subsample is UV bright and has z < 0.4. The 'FUSE-HST' 
subsample has 24 objects, 17 of which come from the Far 
Ultraviolet S pectroscopic Explorer (FUSE) AGN program 
(|Krisd l200lh . This is a heterogeneous, UV-bright sample 
with z < 0.5. The 'RLQ' subsample includes nearly 50 
quasars originally assembled to study orientation; all mem- 
bers of the sample have similar extended radio luminosity 
which is thought to be isotropic. The blazars originally in- 
cluded in this sample are excluded in the SED atlas be- 
cause of their variability due to sync hrotron emission from 
a beamed jet. Using the method of iMassaro et al] (|201ll ) 
and data available from WISE, we verify that none of the 
remaining objects have a large non-thermal contribution to 
emission in the infrared. 

Infrared coverage comes from Spitzer and near- 
infrared coverage is fro m 2 Micron All Sky Survey 

(2MASS; ISkrutskie et al] I2006T ) photometry. Mid-infrared 
spectroscopy was obtained for 46 objects fro m the Spitzer 
Infrared Spectrograph qRS lHouck et al]|2004l ). The infrared 
spectra cover about 3 — 35 fj,m rest frame for the redshifts 
of the sample. Far-infrared photometry at 24, 70, and 160 
/xm from Multiba nd Imaging Photometer on Spitzer (MIPS; 
iRieke et al] 120041 ) is available for 50 objects, all using pho- 
tometry mode. 

AGN variability, which can be a concern, does not ap- 
pear to be a significant issue in these SEDs. The data in 
the infrared through X-ray region of the SED were col- 
lected between 1991 and 2007. Optical-FUV data were taken 
quasi-simultaneously (within weeks) and infrared emission, 
which arises from a size scale on the order of parsecs (e.g., 
iRaban et ail [2009) or less for the hottest dust, will likely 
vary on a longer timescale of months to years. There are 
no strong discontinuities in the infrared emission that might 
indicate variability, but we investigate it in detail. 



WISE data are available for 34 objects in our sample 
that have IRS spectra from Spitzer, allowing us to determine 
infrared variability on the timescale between the Spitzer and 
WISE missions. IRS coverage is such that we are only able 
to make this comparison in the W3 and W4 bands, at 12 and 
22 (jrn respectively. We use the response curv e s and magni- 
tude zero points for WISE from ljarrett et al] (|201ll) to cal- 
culate WISE magnitudes from our Spitzer spectra and com- 
pare these to the observed WISE magnitudes downloaded 
from the WISE All-Sky Data Release. We apply a flux cor- 
rection to fix a known discrepancy that arise s between red 
and b lue calibrator stars, described in detail in I Wright et al.l 
(|2010i ). to the WISE All-Sky Data Release magnitudes of 
+ 17% and -9% in flux to the W3 and W4 bands respec- 
tively. The resulting WISEs P i izer and WISEAii-sky magni- 
tudes are given in Table [1] 

The narrow distribution of differences between WISE 
and Spitzer magnitudes indicates that these objects have 
varied by less than 5 — 10% during the 4 — 5 years between 
the two missions (these Spitzer data were taken in Cycle 2, 
which began in June 2005 and lasted one year, and WISE 
data were taken during the first half of 2010). The obvious 
exception is 4C 11.69, which is highlighted in Table [T] This 
is a core-dominated object with log R = 1.47 that likely 
has a significant synchrotron contribution to emission in the 
infrared that was not classified as a blazar based on optical 
variability by IShang et al] d201lf) or infrared colors using 
the method of IMassaro et alT i 20 1 ll ) . The strong infrared 
variability is clear evidence that this object is a blazar that 
was missed by other identification methods so we exclude 
this object from t he sam ple for the rest of our analysis. 

IShang et al] (|201ll ) applied two important corrections 
to the data while assembling the SEDs. First, at optical-to- 
far-ultraviolet wavelengths the SEDs suffer from Galactic 
extinction from dust. This is removed wit h the empirical 
mean ex tinction law oflCardelli et al] (1 19891) using the dust 
maps of ISchlegel et al.l (| 19981) . SecondT Tshang et al] (|201ll ) 
corrects for host galaxy contamination at near infrared and 
optical wavelengths, although host galaxy contamination to 
the AGN light is not usually large for their UV/optically 
bright q uasars it can be sig nificant even in higher luminosity 
quasars l|Lusso et al,ll201ol ). 

Note that as a result of correcting the 2MASS points 
for host contamination, the shape of our SEDs will be dif- 
ferent than an un corrected SED in t hat region. Host frac- 
tions measured bv lShang et al] (|201ll ) in the J-band had the 
range 0.08 — 0.94, with an average of 0.40. In some objects 
this can be a significant effect, more so at low luminosities. 

We do not include WISE data in the SEDs. For many 
objects we have IRS and MIPS coverage from Spitzer and 
do not gain new coverage by adding WISE, though we do 
find good agreement between the Spitzer and WISE data 
as evidenced by the low variability between the epochs of 
Spitzer and WISE and the WISE data plotted for one ob- 
ject in Fig. [1] WISE data are available for ~ 20 objects 
with bolometric luminosities that do not have Spitzer data. 
While adding these objects would increase the numbers in 
the sample at some wavelengths, it would also increase the 
uncertainty. This would detract from the strength of these 
SEDs, namely the good coverage and detail, particularly of 
the emission features in the infrared that would be wiped 
out by interpolating between WISE points. 
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Table 1. Comparison between WISE and Spitzer 



UDject 


W 3 Spitzer 

(mag) 


\\7 A 

(mag) 


W3 A11 _ Sky 
(mag) 


WTA 

(mag) 


(mag) 


A W7A 

(mag) 


3C 263 


8.28 


6.21 


8.26 


6.12 


-0.02 


-0.09 


3C 334 


8.53 


6.05 


8.54 


5.99 


0.02 


-0.06 


3C 47 


8.11 


6.01 


8.10 


6.01 


-0.00 


0.00 


4C 11.69 c 


7.48 


5.32 


8.71 


6.59 


1.23 


1.27 


4C 34.47 


7.48 


5.38 


7.53 


5.41 


0.05 


0.03 


IRAS F07546+3928 


6.08 


3.86 


6.16 


3.94 


0.08 


0.08 


PG 0052+251 


7.15 


5.10 


7.15 


5.11 


0.00 


0.00 


PG 0844+349 


6.95 


4.94 


6.87 


4.85 


-0.07 


-0.09 


PG 0947+396 


7.55 


5.66 


7.52 


5.57 


-0.02 


-0.10 


PG 0953+414 


7.56 


5.69 


7.38 


5.65 


-0.18 


-0.04 


PG 1001+054 


7.81 


5.97 


7.65 


5.94 


-0.15 


-0.03 


PG 1100+772 


7.61 


5.69 


7.66 


5.61 


0.06 


-0.08 


PG 1114+445 


6.63 


4.53 


6.64 


4.48 


0.01 


-0.05 


PG 1115+407 


7.71 


5.80 


7.63 


5.69 


-0.07 


-0.10 


PG 1116+215 


6.55 


4.69 


6.48 


4.73 


-0.07 


0.05 


PG 1202+281 


7.29 


5.11 


7.44 


5.07 


0.15 


-0.03 


PG 1216+069 


8.13 


6.40 


8.07 


6.45 


-0.07 


0.04 


PG 1226+023 


4.88 


2.86 


4.98 


3.05 


0.10 


0.19 


PG 1309+355 


6.96 


4.69 


7.02 


4.74 


0.05 


0.05 


PG 1322+659 


7.70 


5.73 


7.71 


5.71 


0.01 


-0.02 


PG 1351+640 


5.63 


3.27 


5.65 


3.26 


0.02 


-0.01 


PG 1352+183 


8.11 


6.00 


7.95 


6.05 


-0.16 


0.05 


PG 1402+261 


6.68 


4.77 


6.65 


4.76 


-0.03 


-0.02 


PG 1411+442 


6.36 


4.62 


6.24 


4.59 


-0.12 


-0.03 


PG 1415+451 


7.41 


5.41 


7.43 


5.38 


0.02 


-0.03 


PG 1425+267 


7.85 


5.78 


7.92 


5.74 


0.07 


-0.05 


PG 1427+480 


8.21 


5.80 


8.23 


5.76 


0.02 


-0.04 


PG 1440+356 


6.30 


4.28 


6.31 


4.19 


0.01 


-0.09 


PG 1444+407 


7.33 


5.39 


7.32 


5.35 


-0.00 


-0.04 


PG 1512+370 


8.17 


6.14 


8.24 


6.09 


0.07 


-0.05 


PG 1543+489 


7.04 


4.91 


7.12 


4.82 


0.08 


-0.08 


PG 1545+210 


7.99 


5.94 


7.99 


5.96 


0.00 


0.02 


PG 1626+554 


8.22 


6.75 


8.06 


6.66 


-0.16 


-0.08 


PG 1704+608 


7.08 


4.83 


7.15 


4.83 


0.07 


-0.00 


Mean d 

Standard Deviation 










-0.01 
0.08 


-0.02 
0.06 



a The WISE W3 and W4 bands are at 12 and 22 /im, respectively. 

^ The change in magnitude is taken as the difference in the WISE and Spitzer magnitudes, 
specifically, WISE A11 _ Sky -WISEs p! ( zer . 

c Large values of AW3 and AW4 for this object indicate significant variation of infrared 

emission and a likely synchrotron component. 

™ 4C 11.69 is excluded when computing these statistics. 



2.2 Bolometric luminosities 

We c alculated bolometric luminosities in iRunnoe et al.l 
l|2012l ). Of the total 85 objects in the atlas, only 63 have the 
appropriate wavelength coverage defining a smaller bolomet- 
ric luminosity sample with 23 RQ and 40 RL quasars. The 
RQ quasars are all from either the PGX or FUSE-HST sam- 
ples and are lower redshift (z < 0.5), while the RL quasars 
come primarily from the RLQ sample and more than half 
have higher redshifts [z > 0.5). The RL quasars have an 
average luminosity about 6 times higher than those that are 
RQ. Both RL and RQ quasars span approximately 2 orders 
of magnitude in luminosity. 

Bolometric luminosities were calculated by integrating 
the SEDs from 1 /im to 8 keV under the assumption of 



isotropy. The 1 fj,m lower limit was chosen specifically to 
avoid double counting infrared photons under this assump- 
tion and the 8 keV upper limit was chosen perforce because 
the data do not extend past this limit. Gaps in data cov- 
erage in the near-infrared and extreme-ultraviolet were cov- 
ered by interpolating a log-log power-law spectrum between 
data points on either side. The resulting bolometric lumi- 
nosities have the range log(Lf, i) = 45.1 — 47.3 and are listed 
in Table [2] Uncertainty in the bolometric luminosity has the 
potential to be large, 30 — 40% on average, and is likely sys- 
tematic as well. Since the intrinsic SED in that region is 
unobserved the precise value of the uncertainty is unknown. 
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Figure 1. The rest-frame SED of PG 1202+281. The data are shown by black lines or black points for MIPS and 2MASS. The red 
dotted-dashed line shows the interpolations used to bridge gaps in data. The locations where monochromatic luminosities were measured 
are denoted by dashed gray lines. WISE data points for this object are shown by solid gold triangles with horizon t al bar s indicating the 
frequency range of the WISE filter. The W3 and W4 WISE fluxes have been corrected according to Wright et al. 



2.3 Monochromatic infrared luminosities 

We measured infrared monochromatic luminosities at 1.5, 2, 
3, 7, 12, 15, and 24 microns. These wavelengths are chosen 
to enable comparisons between our bolometric corrections 
and those in the literature and to facilitate the use of our 
corrections for objects observed with WISE or Spitzer at a 
range of redshifts. 

For measuring infrared monochromatic luminosities, we 
defined a smaller infrared sample of 37 objects that have 
data cove rage from IRS a nd M IPS and bolometric luminosi- 
ties from iRunnoe et al.l l|2012t ). In the infrared sample, 23 
objects are RQ and 14 are RL and the RQ objects have a 
factor of 2 higher luminosity on average, though both RL 
and RQ quasars still span approximately 2 orders of magni- 
tude in luminosity. 

We compare our samples in observed parameter space to 
the Sloan Digit al Sky Survey (SPSS) Data Release 7 (DR7) 
quasar catalog l|Schneider et al.l l20T0t ) in Fig. [2] We calcu- 
lated approximate SDSS g magnitude for our objects by con- 
verting the flux at the center of the g-band to a magnitude 
with the SDSS flux zero point. SDSS quasars that have cov- 



erage in WISE and 2MASS are highlighted. In WISE, the 
bolometric luminosity sample covers the magnitude range 
8.4 - 15.5, 7.4 - 14.0, 5.1 - 11.1, and 2.9 - 8.8 at 3.4, 4.6, 
12, and 22 microns, respectively. The bolometric corrections 
derived here will be most accurate for objects that lie in 
the region of parameter space occupied by the sample from 
which they were derived and may be less reliable when used 
on dissimilar objects. 

There were gaps in data coverage in the infrared that 
had to be filled before monochromatic luminosities could be 
measured. Gaps bewteen 2MASS data points, between the 
longest wavelength 2MASS point and the beginning of the 
IRS spectrum, and between the end of the IRS spectrum and 
the next MIPS photometry point were filled by interpolating 
a power-law spectrum between the existing data. 

These interpolations can be seen in Fig. [TJ along with 
the locations selected for making bolometric corrections for 
one representative SED. 

A log-log power-law spectrum may not perfectly de- 
scribe the emission in the regions without data, so there 
is some uncertainty associated with this treatment of the 
SED. We quantify the uncertainty in this step by comput- 
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Figure 2. The region of observed parameter space occupied by 
the 62 objects from the bolometric luminosity sample and the 
SDSS DR7 quasar catalog with the infrared sample distinguished. 
SDSS quasars are small gray points and SDSS quasars with WISE 
and 2MASS coverage are small red points. RL quasars from this 
work are solid blue squares and RQ quasars are open black circles, 
with the infrared sample distinguished by larger symbols. The 
bolometric corrections derived here may lose accuracy for objects 
that occupy another region of parameter space. 

ing monochromatic luminosities again using a WISE com- 
posite (S. Cales, private communication, 2012) scaled to the 
short-wavelength end of the IRS spectrum instead of a log- 
log power-law interpolation. On average, the percent dif- 
ferent between monochromatic luminosities at 3 /im mea- 
sured using the WISE composite and the log-log power-law 
spectrum is 30%, with the luminosities estimated from the 
WISE composite being smaller. The WISE composite is vi- 
sually indistinguishable from a log-log power-law spectrum 
in the region where we make our interpolation, so this uncer- 
tainty is likely an upper limit due to the fact that individual 
objects have different slopes than the WISE composite in 
the gap region. For objects with slopes similar to the WISE 
composite in the gap region, the percent difference between 
monochromatic luminosities at 3 /im using the WISE com- 
posite and the log-log power-law spectrum is on the order 
of a few percent. 

Contamination from the host galaxy can be a con- 
cern at infrared wavelengths. IShanget all (|201lh corrected 
the 2MASS photometry for host contamination. Without 
a prescription for separating the host and AGN emission 
at longer infrared wavelengths, we assumed that the SEDs 
were dominated by the AGN for the corrections at 3 — 24 
fim. This assumption is likely a good one as the host compo- 
nent does not contribute significantly t o the total emission 
until wavelengths lo nger than 24 (jNetzer et alJ 120071 : 
iMullanev etld]l201ll ) 

The monochromatic luminosities, expressed as \L\, 



Table 3. Bolometric ratio statistics 



Wavelength 


Median 


Mean 


ota. ljcv. 


1.5 fim 


12.2 


13.8 


8.9 


2 fim 


10.3 


10.5 


4.3 


3 fim 


9.2 


9.4 


4.2 


7 fim 


9.8 


10.1 


5.2 


12 fim 


9.0 


9.7 


5.7 


15 fim 


9.5 


10.5 


7.2 


24 fim 


10.5 


13.9 


13.1 



were then measured from the SED. The fluxes are well 
known and so uncertainties are small on these measure- 
ments. Monochromatic luminosities at 1.5 fim are measured 
for 62 objects in the bolometric luminosity sample and at 2, 
3, 7, 12, 15, and 24 fim for 37 objects. 



3 DERIVING INFRARED BOLOMETRIC 
CORRECTIONS 

In this section we discuss the derivation of bolometric cor- 
rections cast in several mathematical forms. We determined 
the ratio QraUo = Li SO /\L\ for each object and fit lines to 
log(Li SO ) versus log(AL,\) to determine bolometric correc- 
tions with zero intercepts, Li SO = £ XL\ , and with nonzero 
intercepts, log(L; so ) = A + B log(\L\). We make all three 
corrections at 1.5, 2, 3, 7, 12, 15, and 24 /mi. See Sj4]for a 
discussion about which bolometric corrections are the best 
to use. 



3.1 Ratio bolometric corrections 

The ratio of bolometric to monochromatic luminosity, 
Cratio = Li so /XLx, is the most traditional form of the bolo- 
metric correction. This ratio is listed in Table [2] for each 
object and the distributions are given in Fig. 

The dispersion of infrared bolometric corrections is 
greater than optical/UV bolometric corrections, where the 
bolometric corrections are £ = 4.2 ±0.1 at 1450 A and 
C = 8.1 ± 0.4 at 5100 A, but less than 2-10 keV X-ray bolo- 
metric corrections, where the corrections are £ = 23.31L3.91 
at 2 keV for RL objects and £ = 88.99 ± 30.18 at 2 keV 
for RQ objects. The distributions are fairly similar between 
the infrared wavelengths, except for the few outliers at 24 
/jm. The large 24 fim bolometric corrections tend to result 
from objects that have particularly low 24 fim luminosities. 
Statistics for these distributions are summarized in Table [3] 

We fit a line to log(Li SO ) versus log(\L\) in order to 
determine a bolometric correction with a zero intercept. A 
bolometric correction of this form assumes a single SED 
shape scaled to match a monochromatic luminosity but, be- 
cause quasar SEDs vary in shape, there is no reason to as- 
sume a direct linear relationship with a zero intercept. For 
this reason, we also fit a log-log power law with a nonzero 
intercept to log(L; so ) versus log(\L\). 

Fitting in this paper was done by minimiz ing the chi- 
squared statistic using mpfit (|Markwarddl2009T ) which em- 
ploys the Levenberg-Marquardt least squares method. The 
uncertainties in A La are negligible compared to those in 
Li ao because the spectra are high signal-to-noise and the 
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Table 2. Sample data 

Object Redshift log(/yj so ) Cratio,1.5^tm Cratio,2^tm Cratio,3^tm Cratio,7^im Cratio, 12^m Cratio, 15^im Cratio,24^m 



ergs s 1 


3C 215 


0.4108 


45.77 


11.03 














3C 232 


0.5297 


46.36 


7.21 














3C 254 


0.7363 


46.17 


4.97 


5.36 


4.66 


4.46 


3.71 


3.76 


6.02 


3C 263 


0.6464 


46.81 


10.85 


10.11 


9.20 


10.29 


9.41 


9.76 


10.52 


3C 277.1 


0.3199 


45.61 


14.76 














3C 281 


0.6017 


46.23 


12.81 














3C 334 


0.5553 


46.44 


10.50 


10.32 


9.99 


9.49 


5.79 


5.16 


4.61 


3C 37 


0.6661 


46.18 


8.96 














3C 446 


1.4040 


47.01 


1.64 














3C 47 


0.4250 


45.97 


10.32 


7.28 


4.41 


3.59 


3.32 


3.32 


4.31 


4C 01.04 


0.2634 


45.44 


5.39 














4C 06.69 


1.0002 


47.30 


9.12 














4C 10.06 


0.4075 


46.48 


18.53 














4C 19.44 
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Figure 3. Histograms of bolometric corrections at 1.5, 2, 3, 7, 12. 15, and 24 fim. There are 62 objects total for corrections at 1.5 fim 
and 37 objects for corrections at all other wavelengths. The histogram for 24 fim has much larger dispersion due to several object with 
low 24 fim luminosities. 



monochromatic fluxes are well known whereas the bolomet- 
ric luminosities likely have large uncertainties due to the in- 
terpolation over the unobserved extreme-ultraviolet region 
of the SED. 

Adding a nonzero intercept to the fit is a marginal to 
significant improvement depending on the wavelength, with 
t-ratios given in Table|4] At 1.5 and 24 fim the improvement 
is significant, whereas at 3 — 15 /im the improvement is more 
marginal and at 1.5 fim the nonzero intercept may not be 
needed. 



The fits are displayed in Figs. [3] and [5] and the bolo- 
metric corrections are given in Table [4] Uncertainties in the 
table are 1-sigma uncertainties in the fit coefficients. 

3.2 Radio class 

It is important to justify the use of a single bolometric cor- 
rection for RL and RQ objects, which are known to have 
differences in their SEDs. At infrared, optical, and UV wave- 
lengths the SEDs of RL and RQ quasars are typically very 
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log(AL x ergs s" 1 ) log(AL x ergs s" 1 ) 

Figure 4. log(L(, ;) versus AL A with \ 2 fits f° r 1-5, 2, 3, and 7 /im. The dashed-dotted black lines indicate bolometric corrections with 
zero intercepts, and the solid red lines indicate bolometric corrections with nonzero intercepts. Filled blue squares indicate RL objects 
and op en black circles ind icate RQ objects. 62 objects were fit at 1 .5 (im and 37 w e re fit at all other wavelengths. The dashed gray lines 
are the lEIvis et al ] Gill) bolometric correction at 1.5 /^m and the lRichards et"al bolometric correction at 3 fim. 



Table 4. Infrared Bolometric Corrections 



Wavelength Bolometric correction with zero intercept Bolometric correction with nonzero intercept Sig. of nonzero intercept a 



1.5 ,um L iso = (11.8 ± 0.9) AL A log(L iso ) = (8.98 ± 2.03) + (0.82 ± 0.05) log(AL A ) 4.42 (100.0%) 

2 ^tin L iso = (9.6 ± 0.7) AL A log(L iso ) = (1.85 ± 3.18) + (0.98 ± 0.07) log(AL A ) 0.58 (<44.8%) 

3 fim L iso = (8.4 ± 0.7) AL A log(L iso ) = (4.54 ± 3.42) + (0.92 ± 0.08) log(AL A ) 1.29 (79.8.0%) 
7 ^tm L iso = (8.8 ± 0.8) AL A log(L iso ) = (1.85 ± 4.03) + (0.88 ± 0.09) log(AL A ) 1.57 (88.1%) 

12 ^tm L iso = (8.1 ± 0.9) AL A log(L iso ) = (8.92 ± 4.30) + (0.82 ± 0.10) log(AL A ) 2.08 (95.7%) 

15 nm L iso = (8.5 ± 1.0) AL A log(L iso ) = (10.51 ± 4.39) + (0.79 ± 0.10) log(AL A ) 2.38 (97.8%) 

24 ^tm L iso = (10.1 ± 1.4) AL A log(L iso ) = (15.03 ± 4.77) + (0.69 ± 0.11) log(AL A ) 3.15 (99.6%) 



a The approximate significance of adding a nonzero intercept is given by the t-ratio with the approximate associated probability given in 
parentheses. 
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Figure 5. log(Lj, ;) versus AL^ with \ 2 fits for 12, 15, and 24 jim. The dashed-dotted black lines indicate bolometric corrections with 
zero intercepts, and the solid red lines indicate bolometric corrections with nonzero intercepts. Filled blue squares indicate RL objects 
and open black circles indicate RQ objects. 37 objects were fit at these wavelengths. 



similar, but at X-ray energies they distinguish themselves. 
Because the X-rays are included in the bolometric luminosi- 
ties, this may be a concern. 

The differences between the X-ray emission of RL and 
RQ quasars are well known. For a similar optical flux, 
RL quasars are typica lly three times more X-ray luminous 
|Zamorani et al.l Il98ll ) than RQ quasars but can have an 
X-ray excess compared to RQ quasars as high as a factor of 
3.4 — 10.7 for extr emely radio-loud quasars with R* > 3.5 
(|Miller et al.ll201ll ). 

Despite these differences, in iRunnoe et al.l (|2012l ) we 
found that RL and RQ optical/UV bolometric corrections 
were consistent within the 95% confidence intervals. 

In order to investigate whether separate RL and RQ 
bolometric corrections are required in the infrared, we made 
fits of bolometric luminosity versus monochromatic luminos- 
ity to the RL and RQ objects separately. Because the RL 



objects have higher luminosities on average for the bolomet- 
ric luminosity sample, we isolated radio-loudness by fitting 
only the objects in the region where RL and RQ overlap in 
luminosity. This is not an issue in the infrared sample. 

RL and RQ infrared bolometric corrections are consis- 
tent within the 95% confidence intervals at all wavelengths 
considered in this investigation. Fig.|6]shows fits to bolomet- 
ric versus monochromatic luminosity for RL and RQ objects 
of overlapping luminosity at 1.5 /jm, which has the most re- 
strictive confidence intervals due to the larger sample size. 
We point out, because the slopes for RL and RQ objects 
do appear somewhat different in Fig. |S] that, while they 
are consistent within the 95% confidence intervals, the con- 
fidence intervals are fairly large for such a small sample and 
a chance remains that separate corrections would be appro- 
priate. 

Our result is consistent with the result from 
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Figure 6. Fits to bolometric versus monochromatic luminosity 
at 1.5 fim for RL and RQ objects of overlapping luminosity. RL 
points in the region of overlapping luminosity are blue squares 
and RQ points in the region of overlapping luminosity are black 
circles. Points outside of the region of overlapping luminosity are 
black crosses and are not included in the fits. 95% confidence 
intervals are shaded in blue (solid) for the RL fit and in gray 
(hashed) for the RQ fit. RL and RQ bolometric corrections are 
consistent within the 95% confidence intervals at 1.5 fim. The fit 
itself is a solid blue line for RL objects and a dashed-dotted black 
line for RQ objects. 

iRunnoe et"all (|2012| ) that RL and RQ optical/UV bolomet- 
ric corrections are consistent within the 95% confidence in- 
tervals and we list here the same caveat. The unknown emis- 
sion in the extreme-ultraviolet and soft X-rays was given the 
same treatment in RL and RQ quasars, but if they have in- 
trinsically different emission a different interpolation would 
be required. 

The blazars originally included in this sample have been 
removed and these bolometric corrections are not appro- 
priate for use on blazars. Extreme caution should be used 
for radio-loud quasars expected of being blazars because of 
the higher likelihood of strong contamination from a jet- 
synchrotron component. These bolometric corrections are 
likely not appropriate for such objects. 



4 DISCUSSION 

The infrared bolometric corrections in Table [4] yield bolo- 
metric luminosities that are consistent with bolometric lumi- 
n osities derived from th e optical/UV bolometric corrections 
ofH 

unnoe et al.l (|20 12T ) . The relationship between the bolo- 
metric luminosities estimated from IR and optical/UV cor- 
rections is consistent with a one-to-one relationship within 
the errors. On average, the scatter between bolometric lumi- 
nosities estimated from the infrared corrections with nonzero 



intercepts and bolometric luminosities estimated from the 
optical/UV corrections with nonzero intercepts is 0.21 dex. 
The scatter for bolometric luminosity derived from bolomet- 
ric corrections with zero intercepts is similar. 

Fig. shows a comparison between optical/UV and in- 
frared bolometric corrections. Of the optical/UV bolometric 
corrections, the dispersion is minimized at 1450 A and maxi- 
mized for 5100 A, so bolometric luminosities estimated from 
those corrections with nonzero intercepts are shown versus 
bolometric luminosities derived from the 12 fim correction 
with a nonzero intercept. Comparisons between other wave- 
lengths are similar. 

There are 7 objects in the lShang et al.l (|201ll ) atlas that 
are not included in the bolometric luminosity sample be- 
cause they lack X-ray data, but that do have optical/UV 
and infrared coverage. These objects were not used to de- 
rive any bolometric corrections. For these objects we mea- 
sure monochromatic luminosities and apply the optical/UV 
and infrared bolometric corrections. We find good agreement 
between bolometric luminosities estimated from optical/UV 
and infrared corrections in these objects. 

The bolometric corrections derived here typically yield 
smaller bolometric luminosities than bolometric corrections 
in the literature. For the specific corrections that we discuss, 
this is likely primarily the result of double counting the in- 
frared when integrating the SED to determine bolometric 
luminosity. Emission from the accretion disc is anisotropic, 
so for bolometric luminosity calculated under the assump- 
tion of isotropy, reprocessed infrared photons should not be 
incl uded when integra ting the SED. 

lElvis et alj (1 19941 ) derives the median bolometric cor- 
rection and standard deviation Ltoi = (24.8 ± 8.4) \L\ 
at 1.5 fim. This bolometric correction will overestimate 
bolometric luminosity by a factor of 2 compared to our 
1.5 fim bolometric correction with a nonzero intercept for 
XL\ = 45.0, which is in the middle of the range measured 
for our sample. At the same monochromatic luminosity it 
will overestimate bolometric luminosity by a factor of 2.1 
compared to o ur 1.5 fim b o lomet ric correction with a zero 
intercept. The lElvis et al.l (| 19941 ) bolometric correction is 
not consistent with our 1.5 /im bolometric correction with a 
zero intercept within the g iven errors. 

IRichards et~ai1 £2006) derives a bolometric correction 
and standard deviation Ltoi = (9.12 ± 2.62) \L\ at 3 fim. 
This bolometric correction will overestimate bolometric lu- 
minosity by a factor of 1.1 compared to our 3 /im bolometric 
corr ection with a nonzero or zero intercept for \L\ — 45.0. 
The IRichards et al.1 l|2006l ) bolometric correction is consis- 
tent with our 3 fim bolometric correction with a zero inter- 
cept withintheg^enerrors . 

iHopkins et al.l (|2007l ) derives a luminosity dependent, 
double power law bolometric correction at 15 /im. We mea- 
sure XL\ at 15 /xm for our sample in order to make a di- 
rect comparison. Fig. [5] shows the rati o Cratio = Lj so /\L\ 
at 15 ^m versus log(Z/{, ; ) (similar to iHopkins et al.1 120071 
fig. 1), with their double power-law bolometric correction 
over-plotted. Our data appear consistent with a luminos- 
ity dependent bolometric correction, but the scatter is large 
enough and the range in Lbol small enough that they by no 
means require it. 

The dispersion in bolometric corrections is created by 
variation between individual SEDs, indicative of varying 



© 2012 RAS, MNRAS 000, [TV?? 



Infrared Bolometric Corrections 11 




47.5 



44.5 45 45.5 46 46.5 47 
l°g(L lao ) (12 Mm, ergs s" 1 ) 



47.5 



46.5 



o 



.45.5 



44.5 




44.5 



45 45.5 46 46.5 47 
log(Liso) ( la M m . er S s s_1 ) 



47.5 



Figure 7. Comparison of bolometric luminosities derived from optical/UV bolometric corrections of lRunnoe et al! I <2012h and infrared 
bolometric corrections. The left panel shows bolometric luminosities derived from the 1450 A correction with a nonzero intercept versus 
the 12 [im correction of the same type. The right panel shows instead the 5100 A bolometric correction. In both panels, the dashed red 
line shows where the bolometric luminosities are the same. Open triangles show the good agreement between optical/UV and infrared 
bolometric corrections for objects not included in determining the corrections. 



physical properties between objects. Some of the spread in 
optical/UV bolometric corrections is due to viewing angle, 
this may also be true of infrared bolometric corrections. The 
additional dispersion in infrared bolometric corrections be- 
yond what is seen in the optical/UV may be due to variation 
in torus covering fraction. If true, this indicates that there 
is a small range in torus covering fractions because there 
is only a small excess in dispersion of infrared bolometric 
corrections. 

Because there are many bolometric corrections avail- 
able in the literature, selecting the most appropriate one 
can be confusing. Here we present our recommendations for 
selecting bolometric corrections. 

• Bolometric corrections at different wavelengths are not 
all created equal; there is more dispersion around the mean 
correction at some wavelengths. In some cases, the data do 
not allow a choice between bolometric corrections at differ- 
ent wavelengths, but when they do we recommend that the 
bolometric correction with the least dispersion be selected. 
This means that UV corrections are preferred, then opti- 
cal corrections, infrared corrections (with those at shorter 
wavelengths being preferred), and finally X-ray corrections. 
However, in cases where optical/UV reddening can be iden- 
tified it is better to use an infrared bolometric correction. 

• Bolometric corrections at a given wavelength come in 
different mathematical forms. We recommend bolometric 
corrections with nonzero intercepts because regression anal- 
ysis shows that adding a constant is usually significant. We 
provide bolometric corrections with zero intercepts primar- 
ily for the purpose of comparing to the literature. 

• It is important to apply bolometric corrections to ob- 
jects that are in the same region of parameter space for 
which the correction was derived. Bolometric corrections 
from the bolometric luminosity sample of 62 objects are ap- 



propriate for objects with z = 0.03 — 1.4 and log(-Lf,oi) = 
45.1 — 47.3. Bolometric corrections from the infrared sam- 
ple of 37 objects are appropriate for z = 0.03 — 0.74 and 
log(Lboi) = 45.1 — 47.0. For higher redshift or luminosity 
objects, it m ay be more app r opriat e to use the bolometric 
corrections of [Richards et all |2006f) . 

• We recommend the additional correction to bolometric 
luminosity for the assumption of isotropy. For a randomly 
selected samp le viewed at the angle of 31° suggested by 
iBarthell (| 19891 ) , bolom etric luminosity will be over-e stimated 
by about 33 per c e nt |Nemmen fc Brothertonll2010l . fig. 11.) 
In [R unnoe et al. 

I i|2012l ). we suggested the following cor- 
rection, which should be applied after using infrared, op- 
tical/UV, or X-ray bolometric corrections: 

Lboi = / Liso ~ 0.75 Li SO (1) 

where / is a factor describing the bias of an average view- 
ing angle due to the anisotropic emission from a disc. For 
objects to which the line of sight is known, a theoreti- 
cal c orrection as a function of view ing angle can be found 
from lNemmen fc Brothertonl (2010), fig. 11. Their predicted 
range for this correction to Li SO is approximately 67 to 200 
per cent for a relativistic disc including projection effects 
and limb darkening at viewing angles of zero to ninety de- 
grees. 

Infrared emission is likely to be more isotropic than ac- 
cretion disc emission, but it is still necessary to make this 
correction when using infrared bolometric corrections be- 
cause they were derived from bolometric luminosities that 
included the isotropy assumption. 

• It will usually be necessary to make k-corrections before 
applying bolometric corrections in the infrared. Because of 
the emission features in the infrared, we recommend scaling 
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Figure 8. The ratio ( ra tio = Li BO /\L\ at 15 fim versus log(L(, ;) 
for thi s sample with the lum inosity dependent bolometric correc- 
tion of lHopkins et al. ] 1 120071) with error shaded in gray. The data 
are consistent with a luminosity dependent bolometric correction, 
but due to large dispersion and a small range in bolometric lumi- 
nosity, one is not required by the data. 



the appropriate composite SED from lShang et all (|201lT ) for 
filling large gaps in data rather than interpolating a log-log 
power-law spectrum between available data points. Keep in 
mind when doing this, that the composite SEDs have had 
corrections for host galaxy emission applied to the 2MASS 
points in the near-infrared and to the optical spectrum. This 
will change the shape of the SEDs at these wavelengths com- 
pared to uncorrected SEDs. 



5 SUMMARY 

We have added infrared bolometric corre ctions to the s uite o f 
bolometric corrections derived from the lShang et all (|201lh 
SEDs. We provide corrections with a zero and nonzero in- 
tercept for 1.5, 2, 3, 7, 12, 15, and 24 fim. We find that 
corrections derived for RL and RQ objects are consistent 
within the 95% confidence intervals and, while the data may 
be consistent with a luminosity dependent correction, they 
do not require one. 

The bolometric corrections derived here are self con- 
sistent with the o ptical/UV bolometric corrections of 
iRunnoe et al]ll2012h . The scatter in infrared bolometric cor- 
rections is larger than in optical/UV bolometric corrections, 
but smaller than in X-ray bolometric corrections. Bolomet- 
ric luminosities resulting from these infrared corrections will 
typically be sm aller than those resulting from those in the 
literature ( e.g.. lElvis et all Il994l; IRichards et al.ll200fi ), al- 
though the IRichards et all ([2006) correction is consistent 
within the errors. 

For objects with properties that fall in the region of 



parameter space covered by this sample (log(L{, i) = 45.1 — 
47.0 for 2-24 /an and log(L 6o! ) = 45.1-47.3 for 1.5 fim), we 
recommend the following infrared bolometric corrections: 

• Recommended Correction at 1.5 fim 

log(L iso ) = (8.98 ± 2.03) + (0.82 ± 0.05) log(AL A ) (2) 

• Recommended Correction at 2 /im 

log(L iso ) = (1.85 ± 3.18) + (0.98 ± 0.07) log(AL A ) (3) 

• Recommended Correction at 3 fim 

log(L iso ) = (4.54 ± 3.42) + (0.92 ± 0.08) log(AL A ) (4) 

• Recommended Correction at 7 (im 

log(L iso ) = (6.31 ± 4.03) + (0.88 ± 0.09) log(AL A ) (5) 

• Recommended Correction at 12 /im 

log(L iso ) = (8.92 ± 4.30) + (0.82 ± 0.10) log(AL A ) (6) 

• Recommended Correction at 15 /im 

log(L iso ) = (10.51 ± 4.39) + (0.79 ± 0.10) log(AL A ) (7) 

• Recommended Correction at 24 /im 

log(iiso) = (15.03 ± 4.77) + (0.69 ± 0.11) log(AL A ) (8) 
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